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Abstract 



We present measurements of circular polarization from rotational spectral lines of 
molecular species in Orion KL, most notably ^^CO (J = 2 — )• 1), obtained at the Caltech 
Submillimeter Observatory with the Four-Stokes-Parameter Spectra Line Polarimeter. 
We find levels of polarization of up to 1 to 2% in general, for ^^CO ( J = 2 — )• 1) this 
level is comparable to that of linear polarization also measured for that line. We present 
a physical model based on resonant scattering in an attempt to explain our observa- 
tions. We discuss how slight differences in scattering amplitudes for radiation polarized 
parallel and perpendicular to the ambient magnetic field, responsible for the alignment 
of the scattering molecules, can lead to the observed circular polarization. We also show 
that the effect is proportional to the square of the magnitude of the plane of the sky 
component of the magnetic field, and therefore opens up the possibility of measuring this 
parameter from circular polarization measurements of Zeeman insensitive molecules. 

Subject headings: ISM: clouds — ISM: magnetic fields — polarization — ISM: molecules 
— ISM: individual (Orion KL) 



Introduction 



Magnetic field studies in molecular clouds are most often conducted through the detection of 
polarization signals, either in molecular lines or dust continuum spectra. Except for Zeeman mea- 



Crutcher et al. 


1999: 


Broean & Troland 


2001 



Falgarone et al.ll2008l ). all other types of observations and analyses provide indirect characteriza- 
tions of magnet ic fields. Perhaps th e onl y technique t hat d oes not rely on polarization measure- 
ments is that of lHoude et al.l (j2000al ) and lLi &: Houdd (|20081 ), which is b ased on the compar i son of 



line widths from coex i stent molecular ion and n eutral species (see also iHoude et al.l l2000bl . 12001 
Tillev k. Balsaralboiol . I2OI1I : IPalceta et al.l[201oh . 
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For a c ouple of decades or so measurements of dust continuum polarization, either in emission or 
absorption jHeileslboOol: iNovak et alJ bood: iLi et al.|[2006l : iMatthews et alibood : lOotson et aPboiol : 
Vaillancourt &: Mattliewsll2012l ). have been a particularly efficient of way of conducting systematic 
studies on the interstellar medium and the effects and role of magneti c fields in the processes 



leading to the formation of stars (iGirart et al 



2006 



Attard et al.1 l2009l ) . Although it has been 



known since the work of IChandrasekhar &: Fermil (jl953l ) that polarization maps can be used to 
provide an estimate of the strength of the plane of the sky component of the magnetic field, recent 
developments have rendered possible refinements of this technique and the further characterization 



of magnetized turbulence and its power spectrum (Houdell2004l : iHever et al.ll2008l : iHildebrand et al 



2009; 


Houde et al. 


2009. 


2011: 



Chitsazzadeh et al. 



201^ . 



Linear Polarization of molecular lines through the so-called Goldreicli-Kylafis effect (jGoldreich fc Kylafis 



1981 



, which was first confirm ed observatiqnally in the envelope of an evolved star some 16 years 
after its theoretical prediction (iGlenn et al.lll997l ). is now routinely detected in star- forming regions. 
This effect specifies the conditions leading to an imbalance in the population of magnetic sub-levels 
responsible for the vr- and u-transitions, which are respectively polarized parallel and perpendicular 
to the orientation of the plane of the sky component of the magnetic field. Depending which sub- 
level populations dominate, the detected linear polarization can be oriented in either directions. It is 
i n principle possible to lift this dege neracy through the observations of several molecular transitions 
(jCortes. Crutcher. &: WatsonI |2005| ) . With the development of array heterodyne receivers it is to 
be expected that this technique will also become an important tool for mapping magnetic fields on 
large scales in the interstellar medium. 

The prediction and observations of circular polarization signals, beyond the measurement of 
the Zeeman effect, have proven even more ch allenging. Bu t impr ovements are happening on the 
observational front. Recent observations by iMunoz et al.l (j2012l ) have successfully revealed the 
detection of continuum circular polarization levels of about 1% at wavelengths of 1.3 mm and 860/im 
in Sgr A*. They attribute the presence of such signals to a co nversion of lin e ar pol arization to circular 
polarization (or Faraday conversion). For molecular lines, I Cotton et al.l (120111 ) have reported the 
clear detection of non-Zeeman circular polarization of SiO maser lines in the asymptotic branch giant 
star IK Tau with the VLBA. They found high levels of circular polarization (often exceeding 10%) 
that can also amount to a significant fraction of the linear polarization simultaneously detected. 
They were unable to explain their observations using alternate models for the conversion of linear 



polarization to circular p olarization based on populations imbalance (e.g., IWiebe &: WatsonI Il998l : 
Deguchi fc WatsonI Il985l ). It thus appears that a novel physical model is needed to explain the 



presence of such polarization signals in molecular lines. 

This is what we endeavor to accomplish in this paper, where we present the recent detection 
of circular polarization signals in rotational spectral lines of molecular species in Orion KL, most 
notably ^^CO (7 = 2—7-1) (Sections[2]and[3]). We also introduce a physical model based on resonant 
scattering in an attempt to explain our observations (Section |4j). We follow with a discussion in 
Section [5] and end with a short summary in Section |6l 
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2. Observations 

The observations discussed in this paper were obtained at the Caltech Sub milhmeter Observa- 



tory (CSO) with the Four-Stokes-Parameter Spectral Line Polarimeter (FSPPol: iHezareh Sz Houde 



20101 ). FSPPol allows measurements of linear and circular polarization signals by the insertion of 
half-wave (HWP) and quarter-wave (QWP) plates in the CSO telescope beam before detection using 
the facility's receivers, which are sensitive to a single state of linear polarization. Although FSPPol 
can operate in the 230 GHz, 345 GHz, and 492 GHz bands with the corresponding receivers, the 
observations presented in this paper w ere all conducted in the 2 30 GHz band. More details on the 



functioning of FSPPol can be found in iHezareh &: Houdd (120101 ) 



The circular polarization observations of ^^CO ( J = 2 — )• 1) at 230.5 GHz presented and dis- 
cussed in this paper were obtained on 23 November 2011 under mediocre conditions (r (225 GHz) ~ 
0.2) and on 5 February 2012 under excellent skies (r (225 GHz) 0.03 - 0.04). Both sets of ob- 
servations yielded similar results; the spectrum stemming from the February 2012 observations is 
shown in Figure [TJ Also discussed are subsequent circular polarization measurements aimed at 
the HON ( J = 3 — )• 2) rotational transition at 265.9 GHz, which were realized on 8 and 9 February 
2012 (r (225 GHz) ^ 0.05), and shown in Figure El Finally, linear polarization of ^^CO ( J = 2 -)> 1) 
observations obtained on 12 February 2012 (r (225 GHz) ~ 0.13) are sh own in Figur e [3l T his linear 



polarization spectrum is in good agreement with the previous result of lGirart et al.l (|2004l ) (see the 
two central panels of their Fig. 1). All observations were pointed to the peak position of Orion 
KL at RA(J2000) = 05'^35'"14!5 and Dec (J2000) = -05°22'30'.'4. The telescope efficiency was 
measured to be ~ 60% from scans on Jupiter, while the pointing accuracy was determined to be 
better than approximately 6". 

In order to minimize instrumental effects in the detection of polarization signals due to pointing 
and calibration errors, a conservative observation method was adopted. That is, the necessary 
integrations for the measurement of linear polarization (at four HWP orientations) and circular 
polarization (at two QWP orientations) were kept short at one minute ON-source, and a temperature 
calibration was done before each of them (see IHezareh fc Houdelboioh . The low levels of polarization 



at the peaks of strong lines (i.e., on the order of, or a few times, 0.1% for ^^CO and HCN) found 
in all these spectra are an indication of the benefits in using this method, and can be taken as an 
approximation for the level of instrument polarization present in these observations. 



3. Results 

Our measurement of circular polarization for ^^CO (J = 2 — )• 1) in Orion KL reveals a clear 
detection of polarization levels of up to approximately 1 to 2% across the spectral line, as can be seen 
in Figure [TJ Shown are the Stokes I spectrum, uncorrected for telescope efficiency, and polarization 
levels (symbols with uncertainty, using the scale on the right) in the top panel, while the Stokes V 
spectrum is displayed in the bottom panel. All polarization data satisfy p > 3ap, where p and dp are 
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the polarization level and its uncertainty, respectively. These relatively low polarization levels are 
still significantly higher than the instrumentation polarization expected with FS PPol. For example 



preliminary measurements of the Zeeman effect on CN (A'^ = 2 — ?■ 1) presented in lHezareh fc Houde 



([20101) show that less than 0.2% of Stokes I is expected to leak into Stokes V. This is consistent 
with the polarization level detected at the peak of our ^^CO ( J = 2 — )• 1) Stokes V spectrum (and 
HCN (J = 3 — )• 2); see below), which should be approximately zero in view of the very high optical 
depth expected for this line. Although the probability of getting false Stokes V Zeeman profiles 
(i.e., in the shape of the velocity-derivative of Stokes /) is non-negligible with such observations, 
our spectrum shows no such obvious pattern. 

Nonetheless, steps were taken to ensure that this detection is not spurious. First, the circular 
polarization of ^^CO ( J = 2 — )• 1) was measured twice, once in November 2011 and again in Febru- 
ary 2012, under different sky conditions with similar results. Second, we performed independent 
observations aimed at the HCN ( J = 3 — ?• 2) transition at the same position, which is also known 
to be strong in Orion KL ( Houde et al. 2000a l]; the result from these observations is shown on the 



left panel of Figure [2l As can be seen from the figure, there is no detection of circular polarization 
up to a level of approximately 0.1% in HCN (J = 3 — )• 2), which is an order of magnitude less 
than the ^^CO ( J = 2 — )• 1) detection of Figure [TJ But most interestingly, we see the presence of 
circular polarization at a level of approximately —2% at the peak of a spectral feature located at 
w — 120kms~^, which is mostly dominated by the N^aKc = 12i,i2 — ^ lli,ii transitions of HNCO 
and HN^^CO in the image band at 262.8 GHz; a close-up of this spectrum is shown on the right 
panel of the same figure. The fact that simultaneous measurements on two spectral features (i.e., 
HCN and HNCO/HN^^qq) ^leld none and a detection is strong evidence that our observations 
of ^^CO (J = 2 — )• 1) are not spurious but result from true circular polarization signals present 
in the spectral line. In the next section we present a physical model that seeks to explain these 
observations. 



4. Analysis 

We consider a molecule immersed in a medium that harbors a magnetic field, which provides 
spatial alignment. We further assume the existence of incident radiation close to, or at, a resonant 
frequency of the molecule and in a state of linear polarization at some angle relative to the orien- 
tation of the magnetic field. We can imagine a situation where a population of a given molecular 
species in a background medium emits radiation that is linearly polarized, for example, through the 
Goldreich-Kylafis effect, which is then incident on similar molecules located in the foreground where 



^Although the frequency of that hne is different than the design frequency of the QWP (i.e., ~ 266 GHz vs. 226 
GHz), it can be shown that the only effect caused by the error in the thickness of the QWP at a given frequency is 
a reduction in sensitivity to incoming circular polarization proportional to the cosine of the thickness error. That is, 
no incident linear polarization signal can be converted into circular polarization with the FSPPol set up. 
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the orientation of the magnetic field has changed. We wih investigate the conditions necessary to 
transform hnear polarization i nto circular polarization. Such a situation has been previously con- 
sidered in the literature (e.g., iDeguchi Watsorj|l985f ). but the mechanism for the generation of 
circular polarization presented below is different. 



4.1. Conversion of Linear to Circular Polarization 

We denote this incident background radiation state by |V'o)) which is linearly polarized at an 
angle 6 with the foreground magnetic field, which we assume located in the plane normal to the 
direction of propagation, for simplicity (this assumption will be dropped later on). This radiation 
state can be decomposed as follows 



iV^o) = ao l^ii) +/3o |n_L) , (1) 

where = cos (0) and /3o = sin(^), while |n||) and |nj_) are n-photon states linearly polarized 
parallel and perpendicular to the foreground magnetic field, respectively, propagating toward the 
observer. These states are orthogonal to one another and normalized. The basic idea is to determine 
whether these two states scatter differently off a foreground molecule in such a manner that a small 
relative phase shift (p' is introduced between them. The scattered radiation state would then become 
(up to a global phase term) 



1-0') ~ ao (l + ¥) hll) +/3o|ra±> 

~ aoe^^'|n||> + /3o|nx), (2) 

where it was assumed that the separate phase shifts for both states are much less than unity 
(therefore <C 1 follows for the last equation). It is expected, however, that a very large number of 
scattering N will occur for a given radiation state as it propagates within a molecular cloud. Each 
scattering will contribute a relative phase shift cp' and the final radiation state becomes 

IV') ~aoe^*|n||> + /3o|n^) (3) 

with (j) = N(j)'. 

We can also define n-photon circular polarization states with 

\n±) = ^{\n\\)±i\ni_)). (4) 
It is easy to show that the level of circular polarization (Stokes) v for the final state is 
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V = |Kn+|V)f -IKn-IV')!!' 

= -2ao/3o sin ((^) . (5) 

We can also define the complete state of linear polarization by introducing complementary 
states of linear polarization oriented at ±45° from the magnetic field direction 

|n±45) = ± |n±)) (6) 

and the (Stokes) parameters q and u with 



9 = |Kn|||V)||'-||(n±|V)f 

= «o-/3o (7) 

U = |Kn+45| V')!!^ - ||('^-45| 



= 2ao/3o cos (0) . (8) 
These can be compared with corresponding parameters for the incident radiation state 

90 = al-^l (9) 
no = 2ao/3o (10) 
^^0 = (11) 

to find that, for the chosen system of reference, linear polarization is being converted from uq to 
circular polarization v, and that the total amount of polarization is conserved. 

The previous example considered the case of strictly forward scattering where the radiation 
states before and after scattering were composed of the same jny) and \n^) states. We can generalize 
the analysis by considering a number of initial states propagating in different directions and/or 
containing different numbers of photons (but of the same frequency) such that 

IV-j) = Oij |%,||> + l^j \nj,±) . (12) 

These states could all potentially scatter in the direction of the observer into |ny) and \n_[_). The 
incident radiation state then becomes 
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with ^a'j + /3?^ = 1. The scattered state, ascribed to j = (i.e., |f^o,||) 
\n±)) retains the same form as equation ([2|) (after normalization) 



I nil) and |no 



when 



where 



1^') ~ a'e^'f'' |n|,) + /3'|n^) 



(14) 



a' 



BCao 



1 

D 



Re < 



Z — / -''II ^„ 




(15) 
(16) 



(17) 



B 



C 



D 



Re 



■/3o 



Im 



1 + 



/ /3fc 



5 



y.ll^Um 



-/3o 



(18) 
(19) 

(20) 



and Re {• • •} and Im {• • •} stand for the real and imaginary parts, respectively. In equations (|15p to 
(|20|) (j)'. II and (j)'- ^ are the scattering amplitudes for the different polarization states, while Oj and /3j 
are assumed complex in general (except for oq and /3o, which are chosen to be real). The scattering 
amplitudes (j)'- y and (j)'- ^ will contain geometrical factors that will account for the different incidence 
angles. It can be verified that the earlier forward scattering case is recovered (up to a global phase 
term) when j = A; = is the only possibility, the imaginary components are zero, and </>j|, (/'Y 1- 
The form of equation ([3]) and the others that follow are therefore still adequate for the more general 



case. 



- 8 - 



4.2. Dielectric Susceptibility 



The approach taken in the previous section may appear counter-intuitive with respect to more 
common types of analyses encountered when deahng with molecular polarizability. More precisely, 
it could seem more natural to investigate the relative phase shift between the scattered and incident 
radiations, due to the ii iteraction with a molecule, through calculations of t he iri duced molecular 
dielectric susceptibility (iGrynberg et alJ 120101 : ICohen-Tannoudji et alJ Il988l . 119771 ) . We will show 
here that this approach cannot account for the effect we will discuss in this paper. 

For a detailed and more realistic treatment of radiation/molecule interactions involving two 
molecular energy levels it is often preferable to use the density operator (or density matrix) a to 
perform the analysis. Under this formalism it is found that the j'th component of the mean electric 
dipole moment dj can be evaluated with 



(d,) = TV {ad,} 

= ^Crfearfj.aft, (21) 

a, 6 

where Tr {• • •} denotes the trace and a\,a = (&| ^\o) is a matrix element, with \h) the quantum state 
for level 6, etc. With the knowledge of ^'^j^ can determine the induced molecular polarization 

component Pj = (^dj"^ /V, with V the volume under consideration, and the dielectric susceptibility 
Xj through 



Po=XjEo,j: (22) 

with Eq the incident electric field. 

By performing these calculations for each linear polarization component, one can determine the 
relative phase shift induced between the two scattered electric field components from the difference 
in the corresponding dielectric susceptibilities. For example, if we respectively have x\\ ^iid xi. foi' 
the susceptibilities parallel and perpendicular to the orientation of the magnetic field responsible 
for the alignment of the molecule, then the relative phase shift ip is given by 



^ ^ (X|| - Xa) ^AL, (23) 

with LJ, c, and AL the frequency of radiation, the speed of light, and the propagation path length 
after scattering, respectively. Such calculations, when applied to the case treated in this paper, 
reveal that ip is several orders of magnitude too small to explain our observations. 
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It is important to note that because selection rules for electric dipole transitions specify that 
dj,ab = when a = b, the only components of the density matrix involved in these calculations are 
off-diagonal elements, i.e., aba for b ^ a (see eq. |21)). As we will see in the following section, the 
resonant scattering process we study is of the second order in the electric field and only appears in 
the diagonal elements of the density matrix. More precisely, if an expansion of the density matrix 
in terms of powers of the interaction Hamiltonian Hi is used (iGrynberg et al.ll201Cl ). then it can be 
shown that the second order diagonal terms are 



(2) 



.lVf^(o)_^(o) 



-Ta(t-t') 



to 



{a\ Hi {t') \b) f {b\ Hi {t") \a) e-(-*«+^*")(*'-*")dt" 

Jto 

rt' 

+ {b\ Hi {t') \a) / (a| Hi {t") \b) e-(''^^"+T^")(*'-*")dt" 

Jto 



dt'. 



(24) 



In equation (j24p Ta is the relaxation rate of level a, oj^a ^md are, respectively, the frequency and 
relaxation coefficient for a transition between states \b) and |a), and t — to is the duration of the 
interaction. Since the interaction Hamiltonian involves the electric dipole moment (i.e.. Hi = — d-E), 
it brings a scattering of radiation where the molecule initially in state \a) is momentarily excited to 
(the virtual) state |6) before settling back to |o). As was mentioned earlier, this process cannot be 
captured in calculations involving dielectric susceptibilities described through equations (|2T]) -( l23]l . 



4.3. Resonant Scattering 



We return to the case treated in Section 14.11 of the interaction between an incident radiation 
state and a single molecule, which for simplicity we assume to be linear (like ^^CO). We concentrate 
on two pairs of photon states of the type and |n_L), and njj^ and |n'_j_) for the incident and 
scattered radiation, respectively. These states contain n and n' photons, while in general n 7^ n' . 
It is the n'-photon states njj^ and that are eventually detected by our system to measure 

polarization. 

Since we seek a process that imparts a relative phase shift between such pairs of states, it 
should be clear that the absorption of a photon followed by a spontaneous emission could not lead 
to the desired effect. Such a process would randomize any relative phase difference between the two 
linear polarization states at the emission stage. We must therefore move to a higher (second) order 
mode of interaction. As was mentioned earlier, the best candidate for this is the resonant scattering 
process where an incident photon is absorbed into a virtual, excited state of the molecule and then 
re-emitted into a scattered radiation state. 
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We now denote the initial and final molecule-radiation states as 



\i) = \a)(^\^j) (25) 
I/) = |a'>®|V'), (26) 

with \a) and \a') the initial and final molecular states, respectively at sufficiently long times before 
and after the interaction, while and \^jJ') are given by 



(27) 
(28) 



In view of our earlier discussion we have already included the phase factor e**^' in equation (|28p (see 
eq. |14|). The scattering amplitude result ing from the interaction with a single molecule can then 



be determined through (using MKS units; iGrvnberg et al. 



2010 



Cohen- Tannoudji et al.lll988l ) 



Sif,e 



T 



'nn' 



^ba'i^ba (q1 d • e^/ (6| d-ei\a) 

UJba-^^- ilba 



(29) 



with 00 {(jj') the frequency of the incident (scattered) radiation, uj^a {'^ba') the resonant frequency 
between the initial state \a) (final state |a')) and virtual state "fab is, once again, the relaxation 
coefficient for a transition between the \h) and \a) states (e.g., for a closed system it equals half the 
Einstein spontaneous coefficient, Ai,a/2, when relaxation results only from spontaneous emission). 
As before, n and n' are the number of photons in the initial and final radiation states, while d 
is the molecular electric dipole moment operator and the unit vector associated to the linear 
polarization states with £ = || or _L. The quantities T and are, respectively, the period of 
interaction bet ween the radiation ari d the molecule and the fiducial volume of quantization for the 
radiation field (iGrynberg et al.ll2010l ). From now on, we will assume that in equation (i29l) the final 
and initial linear polarization states are the same, as the difference in the energy levels between states 
of differing polarization (i.e., the Zeeman splitting) will favor similar initial and final polarization 
states (for a sufficiently long interaction period). This also implies that w^a' = ojhai = uj' , and 
\a') = \a). Incidentally, we recognize in equation (I29p the same type of second-order interaction 
term encountered earlier in equation ( |2T 



It follows that because we expect the relative phase shift to be very small, i.e., 



(/.' = Im{%,||-%x}«l, 



(30) 
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we should resist the temptation to ehminate seemingly unimportant differences. More precisely, if 
we write the frequency of the cr-transitions as 



UJ±=U]Q±LOz, (31) 

with ujQ and ujz the vr-transition frequency and Zeeman splitting, respectively, then we should not ap- 
proximate oj± ~ ujQ on the account that loz <^ ojq. We require (/> = Na(j)' ~ 1 for the linear-to-circular 
polarization conversion effect to be measurable when the incident radiation state is interacting with 
a large number Na of molecules. Evidently the size of the volume of interaction is such that the 
relative phase shift (j)' due to the scattering amplitude resulting from the interaction with only one 
molecule, given by equation (j29p . can be extremely small while potentially still sufficient. Taking 
this into account, we rewrite equation (j29p for the tt- and cr-transitions, respectively involving the 
Inii) and \n±) states, with 



if,\ 



-zsm u 



-i sin^ (i) 



T \J nn' 
1? 2eoh 
T V nn' 
l3 2eoh 



(5o| d\\ \a) 



{h± \ d± I a) 



2 (CJQ ± UJz f 

W (wo ± WZ - W - ilb±a) ' 



(32) 
(33) 



where the summation on the virtual states was removed on the account that one state (with 
z = 0, lb) will dominate independently for each transition because of the strong resonance (note 
that oj ^ UJz S> 7feja for the transitions considered here). We have used this notation for the 
virtual states in equations (|32p and (j33p to underline the different types of transitions (i.e., 
vr-transitions bring no change in magnetic quantum number and obey Amj = 0, while fi-transitions 
verify Amj = ±1; hence the notation). Also, the inclination angle of the magnetic field relative to 
the scattering propagation direction (or the line of sight to the observer) is given by t, from which 
d ■ €£ = di sin (i) for the two states of linear polarization. In the numerical calculations presented 
in the next section we will set t = tt/2, effectively setting the magnetic field in the plane of the sky. 

If we now account for the population of molecules with which the radiation interacts, then we 
must also consider the fact that their spectrum (or velocity) will be spread over some normalized 
distribution function h{uj). We can also substitute T = l/c in the same equations, where / is the 
size of the region of interaction. We then have at the frequency oj of the incident and scattered 
photons 



2 / X INayJu {uj) u' (uj) 



sin (l) 



,ba 



{x — Uj) 



-h ix) dx 







d±,ba 





{x + LVz)'^ (x + iOz - ^) {x - OJz)'^ {x - OJZ - ^) 



{x + uz - U)f + -il 



+ 



h (x) dx \ , (34) 
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where Na is the number of molecules in state 



and u 

2 



is the radiation energy density at 
(i = 0,± and i = ||,_L) 



frequency oj. For numerical calculations the terms 
can advantageously be related to the corresponding Einstein spontaneous emission coefficient with 
(jCohen-Tannoudii et al1ll988h 



ba 



H,ba 



(35) 



Notably, for CO we have 



,ba 



\/2 d_\_^ha and — 2Ah^a- It is therefore apparent from 
equation (p¥|) that any relative phase shift would vanish had we approximated loq ± ujz — ojq. 



4-3.1. Circular Polarization of the CO {J = 2 ^ 1) transition in Orion KL 

Let us now provide an estimate of the importance of the effect for the ^^CO ( J = 2 — )• 1) 
transition in Orion KL. The region of interaction I will be constrained by the lifetime of the tran- 
sition, which is det ermined by the relaxation rate 7;,a, or the effective mean free path of a photon. 



Plume et al.l (l2012l ) determined from observations of CO isotopologues that for Orion KL (in the 
Hot Core) the hydrogen number density is riHa = 10^ cm~^ and the temperature Tex = 150 K, 
which together yield a collisional quenching rate for ^^CO of 



721 ~ (<7f) 1.3 X 10 ^ s 



(36) 



This is more than three orders of magnitude greater than = 7 x 10^^ s^^ and we will adopt 
this value for the corresponding relaxation rate, i.e, 75^ = 72 1 (at the specified te mperature the 
momentum-rate transfer coefficient {av) ~ 10~^^ cm^s ^; see Shull &: Draine 1990 ). Accordingly, 
the associated region of interaction equals w c/721 ~ 2.3 x 10^^ cm. The effective mean path of 
a photon will depend on the absorption coefficient a,^ and the resonant scattering coe fficient at 
the frequency oj of the photon. The absorption coefficient (jRybicki &: Lightmanlll979l ) is given by 



Qco (Tex) 



h{uj) 



(37) 



after integrating over the molecular population of the lower state, with 52 = 5 the degeneracy 
of the upper state, Ei is the energy of the lower state {Ei/k ~ 5.5 K, with k the Boltzmanii 



constant), nco 
and Qco {Tex) -- 



10'' cm^ -3 the density of ^^CO (i.e., a relat ive abundance of ap proximately 10 



54.6 its partition function at Tex = 150 K (jPickett et al.lll998l ) 



The resonant scatte ring coefficient i s obt ained through a similar integration of the resonant 
scattering cross-section (iGrynberg et al.l I2OIOI ) over the molecular population of the lower state. 
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which yields 



Qco {Tez 



47r^72i/i(w) , 



(38) 



with the degeneracy of the lower state gi = 3. If we use a Gaussian distribution for the spectral 
profile at its maximum, i.e., h{uj) = 1/ (-v/27rAa;), with Auj = 9.4 x 10'' rads~^ (i.e., a standard 
deviation of approximately 20 kms~^; see below) and oj = 1.4 x 10^^ rads~^ (i.e., 230.5 GHz) we 
find that auj = 8.1 x 10~^^ cm~^ and auj = 4. 7 x 10"^^ cm~^. The effective mean path is then 



determined through ([Rvbicki &: Lightmanlll979l ) 



~ 1.6 X 10^^ cm. (39) 

We therefore find that both path lengths and /^p a-re of similar sizes and little change would 
occur whether we use one or the other. However, it is the case that ^mp ^ ^ and we therefore 
set I = Imp for the size of the region of interaction in equation (I34p . which we rewrite here for 
12C0 (J = 2^1) 



•^^i^ (^) 



(OW — 7^ TTTT^ ^^3^2 Vu{u})u'{uj) 



QCO {Ta 



(x — Cj) 



-h{x)dx 



{x + UJzf jx + ujz - uj) ^ (x - uz)"^ {x-uz - uj) 

(x + UJz - Ojf + ^21 {x-^Z -i^f 



h (x) dx > . 



(40) 



In this equation the volume of the region of interaction was set to ~ . 



It is interesting to note that the integrals in equation (|40p can be combined and transformed 
such that 



0^i"(o;):^u;isin^ (0/ 



mp 



Qco (TLx) 



\J U (bj) U' {uj) I {uj) , 



(41) 



where 



I{lo) = j ^x^{x-uj) 3 (x - w)^ - 7I1 - wl / (x-a;)^ + 7ii 
+ (x - a;) {u? - 3x^) + 7I1 (3x - w) + w| (x + cj) } ^^-^ 



dx, 



(42) 
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with A 



{x + uz -u}) + ill 



{x-ojz- ojf + 7I1 



. We thus find that the effect, or the relative 
to the square of the magnitude of the plane of the sky component 



phase shift 02i , is proportional 
of the magnetic field from the presence of the term w^sin'^ (t) in equation (|4ip . 



For the 10.4-m CSO telescope with an efficiency of approximately 60%, we can convert the 
antenna temperature to energy density according to 



u{uj) = 5.42 X 10 



-22 ^a721 



2.3 X 10 



-35 ' 



erg cm 



(43) 



while our circular and linear polarization spectra shown in Figures [T] and [3] yield u {lo) ~ u' (w) and 
^ '^•^ ^ ™- ii^® wings where polarization is detected. 

Using the parameter values al ready listed before, a s well as i = 7r/2 a nd ujz ~ 1.3 rads~^ 
(i.e., 0.2 Hz, with B ^ I xnG (|Crutcher et al.l Il999l : iHoude et all bood l and ~ -0.269 



(jGordv fc Cooklll984l )). we numerically integrated equations (j4ip and (|42p over a Gaussian profile 
of the form 



1 i/^ ^-"o V 

h(x) = ^= e H Ac. ; (44) 

^/2^Aw 

chosen to approximately match the width of our observed ^^CO ( J = 2 — )• 1) spectra shown in 
Figures [1] and [3j We also used a linear radiation energy density profile u [uj) similar to h but with 
a peak temperature = 1 K. The result is shown in Figure U] where we find that, for this set 
of parameters, the relative phase shift induced by the resonant scattering process is potentially 
significant for ^^CO ( J = 2 — )■ 1) with a maximum value for (j)2i of approximately 30 rad. It is, 
however, important to realize that these calculations are uncertain by a large amount in view of the 
strong dependency of ^g^^ on some parameters. For example, a close look at equation (j4ip reveals 
that (f)2i varies inversely with the fourth power of the gas density tt-Hj and the third power of the 
^^CO abundance; an increase of only a factor of 2 for these parameters would bring the maximum 
value of (1)21 down to approximately 2 and 4 rad, respectively. 

The behavior of (f)2i is different if the gas density decreases to the point where the region of 
interaction is defined by instead of (see eqs. |36| to |39) and the related discussion). We then 
find that the relative phase shift varies inversely with the second power of and is proportional to 
the ^^CO abundance (the same dependencies apply to other molecules with lower abundances than 
^■^CO). Likewise, (l)2i is also strongly dependent on the excitation temperature (mainly through 
the partition function; with the inverse dependency of the molecular abundance) and the magnetic 
field strength (proportional to its second power through the Zeeman splitting). Whatever the case 
it is apparent that the linear to circular polarization conversion effect is strongly affected by even 
modest changes on a range of parameters, and within their established uncertainties. For the 
present case, although there are several combinations that would allow us to match the strength 
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of the polarization conversion effect found in our calculations to the level of circular polarization 
we observe for ^^CO ( J = 2 — )■ 1) in Orion KL, we will simply reduce the magnetic field strength 
from ImG to 0.1 mG (resulting in a Zeeman splitting ojz — 0.02 Hz) and leave all other parameters 
unchanged. Although this value for the magnetic field strength is reasonable in general for star- 
forming regions at such gas densities, there is no guarantee that it is right for this source. The 
corresponding results can be found in Figure [5l where the expected relative phase shift (j)2i for 
^^CO ( J = 2 — )• 1) based on our resonant scattering model (solid curve; left vertical scale), the 
underlying linear polarization radiation profile u (f ) (broken curve), and the corresponding Stokes V 
spectrum (dot-broken curve; right vertical scale) are shown. The Stokes V spectrum is proportional 
to u {v) sin {(t)2i)i and we set 2ajf3j = 1 (see eq. [5]). Its maximum amplitude is now in line with 
that of our observations shown in Figure [TJ We will discuss Figure [5] in more details, especially its 
shape, in Section [5l 

The previous discussion raises the question as to what we should expect observations to reveal in 
cases where (/)2i^ becomes very high. For example, measurements conducted with ^^CO ( J = 2 — )• 1) 
in molecular clouds will commonly probe gases of significantly lower densities than is the case for 
Orion KL. Although changes in other parameters (e.g., the excitation temperature and the magnetic 
field strength) could somewhat offset the effect of a lower gas density, it is likely that under such 
conditions the relative phase shift (/)2i^ will reach significantly higher values (perhaps as much 
as several hundred radians). The fact that the Stokes V spectrum is proportional to sin (</'2i^) 
will bring strong oscillations as a function of the frequency (or velocity) , which will then result in a 
cancellation of any circular polarization across the spectrum (given a finite spectral resolution). This 
cancellation will likely be accentuated if several incident linear polarization modes are resonantly 
scattered into the telescope beam. Finally, as a result of the Stokes u dependency on cos [4>2i^) 
(see eq. [8]), only linear polarization parallel or perpendicular to the magnetic field will be detected 
(i.e.. Stokes g in eq. [7]). 



4-3.2. Circular Polarization of HCN and HNCO in Orion KL 



Evidently, our model applies equally to other linear molecules and should therefore account for 
the lack of circular polarization in our HCN (J = 3—7-2) spectrum of Figure [2l For that molecule, 
its lower abundance, i. e., ~ 10~^, implies thaio ^ = 1.7 x 10^'^ cm in equation (|34p under similar 
conditions as for ^^CO ( Blake et al ] ll987l : ISchilke et al I booih . leaving the corresponding contribution 



to (j) (w) basically unchanged. The reduction in molecular abundance is practically cancelled out by 
an increase of approximately the same factor in the Einstein spontaneous emission (i.e., ^^2^^ — 
8 X 10^'^ s~^) in the numerator of the same equation. However, the Zeeman sensitivity of HCN is 



Because of the larger value of the Einstein spontaneou s emission coefficient for this transition we set 7hcn — 
721 + A§2^ /2 for the relaxation rate I Grvnberg et al. 2O10l l: the value thus obtained for is slightly smaller than 
that for ^^CO ( J = 2 — > f) where the Einstein spontaneous emission coefficient could safely be neglected. 
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approximately three times less than that of ^^CO (i.e., — —0.0962, from lGordv &: Cooklll984 ). 

which implies a corresponding loss of about an order of magnitude for <j) (uj). Taking all these factors 
into account, as well as a decrease of about a factor of four for u (uj), and the increases in uq and the 
partition function (Qhcn (150 K) ~ 210), we find that the effect for HON ( J = 3 — 2) should be 
approximately two orders of magnitude weaker than for ^^CO (J = 2 — )■ 1), again assuming similar 
physical conditions for the two molecular species. This is consistent with our results of Figure 
[2j Moreover, it is likely that the gas density probed by this molecular species is higher given its 
significant critical density (~ 10^ cm~^). An increase in by only a factor of a few would further 
reduce any circular polarization to even lower levels. 

The case of the HNCO/HN^^CO {Nk^k^ = 12i,i2 ^ transitions is more difficult to 

analyze. These molecules are asymmetric tops and will possess more complicated Zeeman spectra 
than linear molecules like CO and HCN (we concentrate on HNCO in what follows). We know, 
however, that the abundance of HNCO (i.e., 10^^ — 10~^) is less than CO by a factor of four to five 



orders of magnitude (iTideswell et al.ll2010l ) and its Einstein spontaneous coefficient is comparable 



to HCN (7 = 3—7-2) with ^H'^CO = 3 x 10^^ s^"*^ for these transitions. However, the Zeeman 
sensitivity is likely to be significantly higher as its electronic ground state possesses electronic spin. 
More precisely, the upper and lower energy le vels involved in the t ransitions discussed here are part 
of triplet states with electronic spin S = 1 (jPickett et al.l 119981). The associate d electronic spin 
contribution to the Lande factor can be approximated with (jGordv &: Cooklll984l ) 



^HNCO ^ 



J {J + 1) + S {S + 1) - N {N + 1) 
J(J + 1) 



(45) 



which for the strongest of these lines, i.e., when AJ = +1 (jPickett et al.lll998l ). covers a range of 
—0.182 to 0.167. The Zeeman splitting between corresponding tt- and cr-lines can be approximated 
by multiplying these Lande factors by the Bohr magneton. Taking into account the fact that the 
nuclear magneton is used instead of the Bohr magneton for similar calculations with ^^CO (and 
HCN), we expect the Zeeman splitting of the strongest HNCO {NxaKc — ^"^1,12 Hi, 11) lines 
to be roughly a thousand times larger than for ^^CO. This would lead to an approximately six 
orders of magnitude increase in the value of in equation (I34p . Combining these changes (i.e., 
in molecular abundance, Einstein spontaneous coefficient, and Zeeman sensitivity) with those for 
the partition function (Qhnco (150 K) ~ 2800, about 50 times that of CO) and the radiation 
energy density (a reduction of approximately 40 for similar levels of linear polarization), we would 
expect the linear-to-circular polarization conversion effect to be at approximately as strong for 
HNCO (NxaKc = 12i,i2 111, 11) as for ^^CO ( J = 2 ^ 1). Although there is a significant level 
of uncertainty in the previous calculations, this result is also consistent with our observations (Fig. 
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5. Discussion 

The analysis presented in tiie previous section established that the relative phase shift induced 
by the resonant scattering process between linear polarization components parallel and perpendic- 
ular to the plane of the sky component of the magnetic field can account for the levels of circular 
polarization detected in the spectra presented in Figures [T] and [21 It may be surprising that a 
Zeeman splitting on the order of 0.1 Hz or less could be responsible for such a significant effect 
across a spectral line that is on the order of 10 MHz wide and centered at 230.5 GHz. But in the 
case of ^^CO ( J = 2 — t- 1) the conversion of linear to circular polarization is likely to be efficient in 
Orion KL, although we again emphasize that small changes in some of the main parameters (e.g., 
or nco) can significantly affect our results. There is one aspect, however, that requires 
further discussion: the expected profile of the circular polarization Stokes V spectrum. 

As was mentioned earlier, Figure [5] also shows the Stokes V profile that results from the 
numerical calculations of (f>2i' discussed in Section 14.3.11 (the dot-broken curve; using the vertical 
scale in the right side of the graph). We approximated this spectrum with the following function 
(in units of Kelvin) 

V (cj) = e"K^) sin [4>^^ (u)] , (46) 

i.e., we set 2aj(3j ~ 1 (see eqs. [5] and |27| ) and, as is apparent in equation (|46p . used a Gaussian 
profile for u {oj) with a peak antenna temperature of 1 K for the radiation associated with the |n||) 
and \n±) incident states. We also converted the abscissa to a velocity scale to ease the comparison 
with the observed spectrum of Figure [TJ The most obvious discrepancy between the two spectra is 
the fact that our calculations yield an antisymmetric profile, while the observations do not. In fact, 
the calculated profile is not unlike the typical Stokes V spectrum expected from Zeeman sensitive 
molecular species and transitions. A detailed comparison between calculations and observations is 
complicated by the fact that our model contains several approximations. For example, the Gaussian 
profile chosen for our calculations does not perfectly match that of the measured linear polarization 
spectrum of Figure [3l and the core of the observed polarization line profiles (linear and circular) are 
very likely dominated by instrumental polarization, which our model does not consider. But it is 
important to realize that several factors can affect the shape of the calculated Stokes V spectrum: 

1. The perfect antisymmetry seen in the results of Figure [5] stems from the symmetry of the 
underlying Gaussian used for the calculations. This is made clearer in Figure [H] where another 
calculation for the same parameters used for Figure [5] is shown, but with the Gaussian profile 
of equations (|44p and (j46p replaced with a slightly uneven line shape. The result is a circular 
polarization spectrum that is markedly broader on one side (where f ^ 15 kms~^) than the 
other. The same would happen for a typical Stokes V spectrum from a Zeeman sensitive 
transition for a slightly uneven Stokes / profile, since they are linked through a derivative. 
But it is important to note that, in our case, it is the line profile of the incident (background) 
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linear polarization radiation that is in question, which may be different than the Stokes I 
spectrum. 

2. In the calculations of the integral contained in equation (p^ (or |42| ) we assumed that a 
photon at a frequency oo will scatter off all molecules at any other frequencies covered under 
the profile h{oj). We must realize that this is not likely to be the case since the size of the 
molecule-radiation region of interaction is relatively small (i.e., ~ 10^^ cm or ~ 10~^ pc; see 
eq. [39]) whereas molecules belonging to different velocity ranges are likely to be separated by 
larger distances within Orion KL (the ^^CO ( J = 2 — )• 1) spectrum covers a velocity range of 
~ ±50 kms~^). This could significantly affect the shape of the resulting Stokes V spectrum. 

3. The numerical calculations that produced Figure [5] considered a single incident radiation mode. 
But we know from the related discussion in Section [4.11 (see eqs. |13| to |20) ) that even at a 
single frequency several incident modes, coming from different orientations with potentially 
differing linear polarization states, take part in the resonant scattering process. This brings an 
averaging process that will certainly affect the line shape of the circular polarization profile. 

4. There also exists an averaging process due to the size of the telescope beam. At the distance 
of Orion KL (i.e., approximately 450 pc) our telescope beam of ~ 32" covers a region of 
approximately 0.07 pc or 10^^ cm, which is several orders of magnitude larger than the region of 
interaction. Our observations therefore contain contributions from several radiation-molecule 
interaction regions, where key parameters that affect the Stokes V line profile are likely to 
vary. 

Although these considerations reveal the complexity of the problem and that we perhaps should 
not expect our observations to closely match the profile resulting from our simplified model and 
presented in Figure [5l it is nevertheless possible to find relatively simple conditions that would 
allow us to calculate Stokes V profiles that are consistent with our circular polarization detection 
in 12C0 (J = 2^ 1). 

To do so we consider incident radiation modes for which (v) varies with velocity (frequency) 
and possesses some distribution about the orientation of the foreground magnetic field on the plane 
of the sky. Since the linear to circular conversion of polarization will only occur when 6 (v) ^ (and 
7r/2) it is possible that most linear polarization signals at frequencies satisfying this condition will 
be efficiently converted to circular polarization leaving only a well defined orientation angle (i.e., 
9 = 0) for the outgoing linear polarization radiation. Our linear polarization spectrum of Figure [3] 
indeed shows a well defined polarization angle across the spectral line. If the distribution of 9 (v) is 
not uniform about ^ = and, for example, is positive in some regions and negative in others, then 
the shape of the circular polarization spectrum would clearly depart from the asymmetric profile 
displayed in Figure [5] (or even [6]) . This is exemplified in Figure [7] where have computed the Stokes 
V spectrum for conditions similar to those used for Figure [5l with the exception that the angle 
9 (v) between the incident linear polarization and the magnetic field orientation varies linearly with 
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velocity at a rate of dO/dv = —1 deg/(kms~^) (with = at = 0). The resuhing Stokes V 
spectrum fohows from equations jS]) and (j4ip (and eq. |46|). It is then found that the presence of 
the factor Uj = sin [6) in this equation results in a symmetric profile about f = that is consistent 
with our circular polarization measurements on ^^CO (J = 2 — )■ 1). As stated above, the core of 
this spectral line is very likely dominated by instrumental polarization and the true polarization 
level there should be zero or close to zero, as in our Figure [71 Although this does not ensure that our 
model can perfectly account for our observations across the whole spectrum, the results of Figure [7] 
are consistent with our measurements away from the core of the spectral line where signals in the 
Stokes V profile have the same sign. 

We also note that the levels of circular polarization predicted by our model are likely to be 
consistent with the IK Tau S i O v = 1, f = 2, ( J = 1 — > 0) (at 43.1 GHz and 42.8 GHz, respectively) 
observations of lCotton et al.l ()201ll ) . Although a more precise analysis would be necessary to ensure 
this (we are not, for example, accounting for possible maser saturation), we can make an approximate 
assessment of the strength of the effect for SiO relative to CO using our previous calculations. If 
we take the SiO v = 2, (J = 1 — )• 0), 42.8 GHz transition as an e xample, we know that although 
the abundance is one to two orders of magnitude less than ^^CO ( Decin et al.ll2012l ). its Einstein 
spontaneous coefficient is about an order of r nagnitude greater at A :^ 4x 10~^s~^, while its Zeeman 
sensitivity is only a factor of two or so lower (jPavis fc Muenteil (|l974} ) find gf^ ~ -0.154 to -0.155 
for v = 0, 1, and 2). If we add to this the fact that the frequency of SiO v = 2, ( J = 1 — )• 0) is 
less than five times that of ^^CO ( J = 2 — )• 1) and set the spectral line width to match their 
observations (on the order of 1 km s^^), then we find that the polarization conversion effect is likely 
to be important for SiO for a significant range of gas densities and magnetic field strengths. A 
similar result is expected for the SiO f = 1, (J = 1 — >• 0), 43.1 GHz transition. These calculations 
could then resolve the known problem of high circular polarization levels found in maser transitions 
that require un reasonably lar ge magnetic field strengths when interpreted within the context of the 
Zeeman effect (jWatsonll2009l ). Furthermore, our model is also consistent with other observations of 
SiO maser lines in evolved stars that showed a correlation between measured levels of linear and 
circular polarizations, no correlation between levels of circular polarization an d Stokes /, as well as 
a case where circular polarization is more important than linear polarization (jHerpin et al.ll2006l ). 



We stress, however, that although some SiO maser Stokes V line profiles presented in lCotton et al 



(|201ll ) match well the asymmetric profile resulting from our simplest calculations shown in Figures 
[5] and [6] (i.e., when the angle 6 between the incident linear polarization and the magnetic field ori- 
entation is constant across the spectrum), many of their spectra show (quasi-) symmetric Stokes V 
profiles that are different from anything we presented. More precisely, if the ^^CO (7 = 2—7-1) line 
in Orion KL is expected to yield little to no polarization at or near the center of the line, the same 
ca nnot apply to maser lines in view of the nature of the stimulated emission process. That is, several 
of lCotton et al.l (120111 ) Stokes V spectra show a maximum (in the absolute sense) near the center of 
the line, not a near-zero value. It follows that the symmetric Stok es V spectrum of Figure [3 cannot, 
by itself, explain the aforementioned results of lCotton et al.l (|201ll ). It is still possible, however, that 
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the presence of several maser spots within the telescope beam could lead to such (quasi-) symmetric 
line profiles under the assumption that the underlying Stokes V spectrum of a given maser is as 
shown in Figure [7] and the masers' systemic velocities are spread across the extent of the observed 
spectral line. The combination of these relatively-shifted Stokes V spectra (i.e., as in F ig. [7]) could 



thus yield something akin to the (quasi-) symmetric line profile of ICotton et al.l (j201ll ). However, 
it is perhaps safer at this point to await the results of an analysis similar to the one presented in 
this paper, but specifically tailored to maser emission, to find out if our resonant scattering model 
can account for all aspects of SiO maser circular polarization spectra with a minimum number of 
assumptions. This we intend to attempt in a future publication. 

Our results also have implications for Zeeman measurements in non masing environments (e.g., 
in molecular clouds). For example, it is clear from our discussion of Section [4.3.21 that HNCO is a 
Zeeman sensitive molecule (in contrast to CO and HCN). However, it is likely that our observations 
of Figure [2] would be dismissed as due to instrumental artifacts when analyzed within the context of 
the Zeeman effect (e.g., leakage from Stokes / to Stokes V). But when studied in conjunction with 
the HCN result presented in the same spectrum and that for ^^CO ( J = 2 — )• 1) shown in Figure [H 
another interpretation is warranted. 



6. Conclusion 



We presented measurements of circular polarization from rotational spectral lines of molecular 
species in Orion KL obtained at the Caltech Submillimeter Observatory with the Four-Stokes- 
Parameter Spectra Line Polarimeter. We measured levels of polarization of up to 1 to 2% for 
i2co (J = 2 ^ 1) and HNCO/HN^^CO {Nk^k, = 12i,i2 ^ Hi.ii), while none was detected for 
HCN ( J = 3 — 7- 2) . We further presented a physical model based on resonant scattering in an 
attempt to explain our observations, through the conversion of linear polarization to circular polar- 
ization. We also showed that this effect is proportional to the square of the magnitude of the plane 
of the sky component of the magnetic field, and therefore opens up the possibility of measuring this 
parameter from circular polarization measurements of Zeeman insensitive molecules. We intend to 
study this in an upcoming paper. 
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Fig. 1. — Circular Polarization spectrum of the ^^CO (.7 = 2—7' 1) observations made at the peak 
position of Orion KL (RA (J2000) = 05''35™14!5, Dec (J2000) = -05°22'30'.'4) on 5 February 2012 
at the CSO with FSPPol. Top: Stokes / spectrum, uncorrected for telescope efficiency, and circular 
polarization levels (symbols with uncertainty, using the scale on the right). All polarization data 
satisfy p > Sap, where p and ap are the polarization level and its uncertainty, respectively. Bottom: 
the Stokes V spectrum, also uncorrected for telescope efficiency. The frequency resolution of the 
Stokes / spectrum is 61 kHz (0.08 kms~^), while the Stokes V spectrum was smoothed by a factor 
of 20. 
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Fig. 2.— Le/i; Same as Figure [2 but for HCN (J = 3 ^ 2), obtained on 8 and 9 February 2012. 
The spectral feature located at ~ —120 kms~^ is a blend of lines from a few molecular species, 
most notably HNCO and HN^^CO in the {NxaKc — ^'^1,12 Hi, 11) transitions; a close-up of this 
spectrum is shown on the right panel. The frequency resolution of the Stokes / spectra is 61 kHz 
(0.07 kms~^), while the Stokes V spectra were smoothed by a factor of 40. 
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Fig. 3. — Linear Polarization spectrum of the ^^CO ( J = 2 ^ 1) observations made at the peak 
position of Orion KL (RA (J2000) = 05^35'"14!5, Dec (J2000) = -05°22'3(X'4) on 12 February 2012 
at the CSO with FSPPol. Top: Stokes / spectrum, uncorrected for telescope efficiency, and linear 
polarization levels (symbols with uncertainty, using the scale on the right). All polarization data 
satisfy p > 3ap, where p and dp are the polarization level and its uncertainty, respectively, and are 
corrected for positive bias in the polarized flux. Middle: Polarization angle from north, increasing 
eastwards. Bottom: Polarized flux pi, also uncorrected for telescope efficiency, but corrected for 
positive bias due to noise. The frequency resolution of the Stokes / spectrum is 61 kHz (0.08 kms""^), 
while the polarized spectrum was smoothed by a factor of 20. 
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Fig. 4. — Numerical calculations of the expected relative phase shift (t>2p (solid curve) for 
^^CO (J = 2 — > 1) in Orion KL based on our resonant scattering model; using the left vertical 
scale. The magnetic field strength was set to 1 mG. The underlying linear polarization radiation 
profile u (v) is shown with the broken curve, using the vertical scale on the right. 
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Fig. 5. — Numerical calculations of the expected relative phase shift 021^ (solid curve) for 
"^^CO ( J = 2 — > 1) in Orion KL based on our resonant scattering model; using the left vertical 
scale. The underlying linear polarization radiation profile u (v) is shown with the broken curve, has 
been normalized in the figure but has a peak antenna temperature of 1 K. The corresponding Stokes 
V spectrum is also shown (dot-broken curve); using the vertical scale on the right. This spectrum is 
proportional to u (v) sin (</>2i^). For this simulation, and the other that follows, the magnetic field 
strength was set to 0.1 mG. 
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Fig. 6. — Numerical calculations of the expected relative phase shift 02]^ (solid curve) for 
^^CO (J = 2 — 7- 1) in Orion KL (using the left vertical scale) for a slightly uneven underlying 
(normalized) linear polarization profile u{v) of 1 K peak antenna temperature (broken curve). The 
corresponding Stokes V spectrum is also shown (dot-broken curve); using the vertical scale on the 
right. Note the change in the resulting line profile when compared with Figure [5j The magnetic 
field strength was set to 0.1 mG. 
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Fig. 7. — Numerical calculations of the expected relative phase shift (/)2i (solid curve) for 
^^CO (J = 2 — )• 1) in Orion KL (using the left vertical scale) for a case where the angle 9 between 
the incident linear polarization and the magnetic field orientation varies linearly with frequency 
(or velocity) at a rate of dO/dv = —1 deg/(kms~^) (with = at = 0). The underlying lin- 
ear polarization radiation profile u{v), shown with the broken curve, has been normalized in the 
figure but has a peak antenna temperature of 1 K. The corresponding Stokes V spectrum is also 
shown (dot-broken curve); using the vertical scale on the right. This spectrum is proportional to 
— 2 sin {9) cos {9) u {v) sin (021^) (see eq. [5]). Note that the Stokes V spectrum now has a symmetric 
profile about u = 0. The magnetic field strength was set to 0.1 mG. 



